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Abstract 

We show that, contrarily to recent claims, data from the Wilson (unimproved) fermionic action 
at three different values demonstrate the presence of a large Goldstone boson contribution in the 
quark pseudoscalar vertex, quantitatively close to our previous estimate based on the SW action with 
csw = 1.769. We show that discretisation errors on Zg IOM /Zp IOM seem to be much smaller than 
the Goldstone pole contribution over a very large range of momenta. The subtraction of this non 
perturbative contribution leads to numbers close to one-loop BPT. 

LPT-Orsay 00-122 (november 2000) 



1 Introduction 

In a recent paper [Q, using dataf] from the QCDSF collaboration at (3 — 6.0 and csw — 1-769 for three 
different values of k we have shown that the quark pseudoscalar vertex contains an unexpectedly 
large contribution from the Goldstone boson pole, and that this contribution accounts for a third of 
Zp at 2 GeV. Qualitative indication of this large Goldstone contribution had been noted previously 
by the QCDSF group and by the Rome group A large quantitative estimate has been found 
independently by JLQCD Jj| with staggered fermions. 

Giusti and Vladikas |^| have recently presented a criticism of our paper; they have reexamined this 
problem using Wilson data at two values of /3 = 6.2, 6.4, and have come to the conclusion that the 
Goldstone boson term would be below the level of discretisation errors "around p = 1/a", and therefore 
not significant. However, they do not compare their data with our result — given at 2 GeV (ap ~ 1 at 
(3 = 6.0)— but rather consider higher momenta p — 3.3 or 4.6 GeV (sin 2 (ap) = 0.8 in their Fig. 4 and 
Table 1), where of course the Goldstone is much smaller^. We present here an analysis based on a set 
of previously published data of the same origin [[7), which shows that, contrary to their objections, the 
Goldstone boson contribution in Wilson data is in fact completely compatible with our previous estimate, 
and much above discretisation errors at 2 GeV and probably at notably higher momentaQ. 

The present study in fact improves our determination of the Goldstone pole. First of all, the set of 
Wilson data analysed here is obtained at larger values of /?, up to f3 = 6.4, where discretisation errors 
should become really small at moderate p ~ 2 GeV (a 2 p 2 ~ 0.25 at j3 = 6.4), and perhaps better than 
with the previous (3 = 6.0 data Q with ALPHA SW improved action. Secondly, we are now in a position 
to give a reliable estimate of the discretisation errors by considering the evolution of the parameters with 

1 Unite Mixte de Recherche - CNRS - UMR 8627 

2 To calculate Z^, we have also used the propagator data at f) = 6.0 kindly communicated by the Rome group. 

3 One should not speak of the magnitude of the Goldstone at p = 1/a independently of the value of /3, since this magnitude 
denends on p, and p = 1/a depends of course on /3. We have never spoken ourselves of such a thing, as stated in the abstract 
of M|, but of its magnitude at 2 GeV (J3 = 6.0). 

The present study could presumably be improved using the raw data of M], which were unfortunately not available to 
us up to now. We now hope to improve our results in the future with the help of the authors of Jg). 



/3. Thirdly, thanks to the (3 = 6.2,6.4 data, we can improve the large momentum tail of our analysis of 
power corrections. Finally, we shall also improve our previous work by the inclusion of statistical errors. 

Having shown that the Wilson data confirm rather beautifully our first estimate, exhibiting a remark- 
able stability of the effect with increasing 0, we shall conclude by a critical analysis of the procedures of 
H which lead to erroneous conclusions. 



2 Previous results on the Goldstone pole in the pseudoscalar ver- 
tex 

2.1 Our previous results 

The theoretical expectation from the continuum is that a pole in 1 jm q must be present in the pseudoscalar 
quark vertex at q = 0, as a consequence of the existence of the Goldstone boson. In Jl[], we analysed the 
lattice data kindly communicated by QCDSF collaboration [|| for the PS vertex at (3 — 6.0 with SW 
action at csw — 1-769, at several k, combined with propagator data from the Rome group, at the same 
(3 with the same action. We have obtained, in the MOM renormalisation scheme, a Goldstone-like fit of 
(Zp l ^ MOM = Tp/Z^p as function of k. Namely, we have shown that, at ap = 1 and f3 = 6.0, i.e. around 
1.9 - 2 GeV: 

Z P J (2 GeV) = 1.88+—. (1) 
am q 

The first term on the r.h.s. is the (/?-dependent) short-distance contribution. The Goldstone pole cor- 
responds to the second termF], i.e. a pole in m q at m q — 0. We have also checked that, as function 
of p 2 , one has the expected behavior: the short distance term is compatible with a logarithmic depen- 
dence ~ [a s (p 2 )] 4 ^ n and the Goldstone term has a 1/p 2 decrease. Converting to physical units, with 
a -1 = 1.9 GeV, one obtains: 

n,-\, 2n „-i, , , \ 0.158 GeV 3 

Z'p \P ) = Z'p (short distance) H ~ (2) 

m q p 1 



Of course, the effect of the uncertainty due to the error on a -1 could be relevant for the Goldstone 
contribution since it is cx a~ 3 . Despite this fact, and despite the presence of other uncertainties, it 
seems difficult to escape the conclusion that the magnitude of the Goldstone term is large at the smallest 
quark mass (around m q = 50 MeV) and at 2 GeV: 30% of the total Zp 1 — 2.7 at 1.9 GeV, although it is 
decreasing rapidly with increasing p 2 . The result can be translated into an estimate of the Georgi-Politzer 
mass at 1.9 GeV in the chiral limit: rap = 34 MeV. 



2.2 Related findings of JLQCD and ALPHA 

Our evaluation (|j],|j) is quantitatively supported by the remarkable JLQCD results on the pseudoscalar 
vertex and the mass operator with staggered fermions Jj|. These results are important as they benefit 
from two advantages: they go down to very small quark masses (about 20 MeV), and they have, in 
principle, small discretisation errors. The phenomenon appears very stable with respect to (3 as they 
considered (3 = 6.0, 6.2 and 6.4. 

The above estimate of the Goldstone term is also supported by the estimate of the ALPHA group 
for the short distance Zp Q, which must be considered as very solid, since they work at ultra-short 
distances, and since their discretisation errors are very well controlled. When their short-distance result 
is converted into the MOM scheme^] and evolved perturbatively (at 3 loops ) down to 2 GeV, one obtains 
Zp 1 (2 GeV) = 1.8. This result is close to the first term of Eq.(|l|), and quite different from the total 
Zp X {2 GeV) ~ 2.5 — 2.7: the difference must be filled by the Goldstone boson pole, unless there be 
incredibly large discretisation errors in the total Zp 1 . The latter is very unlikely in view of the following 
discussion of Wilson data. 

6 Recently, on investigating the quark propagator and the Ward identity relating the PS and the propagator jfj], we have 
improved the precision of_the determination of Zp, and obtained similar numbers. However, we shall stick here to our first 
determination, to which [rj is referring. 

6 The initial idea of this conversion is due to Vittorio Lubicz. 



2 



3 The Goldstone pole in Wilson data 



The most valuable part of § is the introduction of the ratio (Z P / Z S ) RI /MOM = Zf OM /Zf OM 
and, on the other hand, of some interesting Ward-Takahashi identities. Let us emphasize indeed that 
r?MOM jgMOM - g a gcaie-dependent quantity, in contrast to (Z P /Zs) WI , but with a p 2 dependence due 
only to power corrections . This gives it an important advantage over Z P which necessarily contains 
a purely perturbative contribution with logarithmic behaviour, complicating the determination of the 
power term. 

We therefore discuss Zp IOM /Z¥ as the best probe of the Goldstone pole, which should be seen as 
a l/m q p 2 term in the inverse, Z¥ /Zp . Moreover, we shall show later that, according to equation 
(19) of §, a p 2 change in Z^ IOM /Z¥ OM can only be due to the presence of a Goldstone l/m q pole. 
The physical source of any departure from the p 2 — > oo asymptotic value must then be a Goldstone 
contribution. That it is present is recognized in H; we differ on the estimate of its magnitude. 

The first question to be answered is whether the effect of the Goldstone pole has the large magnitude 
that we have estimated, or whether it is sub-dominant with respect to discretisation errors already at p = 
2 GeV, as claimed in || . This can be answered only by considering the behaviour of z¥ OM jZp OM around 
2 GeV, or, in a scale independent manner, by comparing the coefficient of the power corrections to the 
one we have given in Eq. (||). 

The second point concerns the estimate of the discretisation error itself: its magnitude can be esti- 
mated by examining the stability of the result as a function of with fixed physical parameters. 



3.1 Methodological considerations 

3.1.1 Use of physical units and comparison of different actions 

We have to compare different sets of data, with different actions and different /?'s. It is a delicate task, 
especially for scale-dependent quantities. Since the Goldstone pole residue is a physical effect — although 
perhaps gauge dependent — seen in the renormalised pseudoscalar (PS) vertex, a minimum requirement 
is to compare the results at identical momenta for the same quark mass, not at identical ap if p varies, 
as done in ref. ||, e.g. when making statements about the magnitude of the Goldstone as compared to 
ours, at "p — 1/a". Our Fig. 1 below illustrates the effect of comparing data in terms of ap instead of 
p: the data at various /3's, which show large discrepancies in terms of ap (Fig. 1 a), almost superpose in 
terms of physical p (Fig. 1 b)f\. 
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Figure 1: (a) the data from @ for zf OM /Z^ OM for three values of p, as function of the lattice a 2 p 2 
and (b) our fits to these data represented in physical units. 

7 This is also quite visible in Figs. 3 to 5 of M. 
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The comparison of data from the same action at different (3 requires only the ratio of the lattice units 
a, which are rather well determined. When we consider the absolute magnitude of the Goldstone term, 
the uncertainty becomes larger since the strength (i.e. the coefficient of \/(m q p 2 ) ) is proportional to 
a" 3 . 

Furthermore, a dependence on the action, and therefrom an additional cutoff dependence, is to be 
expected even on the the finite Z's as shown by lattice perturbation theory. In particular, the Wilson 
term and the clover term induce contributions to finite Z's of the form Cg 2 — C6.0//3, with a coefficient C 
dependent on the action. These contributions are quite sizable — even for large (3, close to the continuum. 
On the other hand, one may expect that the non-perturbative Goldstone part is independent of the action, 
since it is a long distance effect. Both these expectations are confirmed in the present analysis. 



3.1.2 Discretisation errors 

At a given /?, one can appreciate discretisation errors, as done in ||, by observing the discrepancies 
between various quantities which should have been equal, for instance between various estimates of 
(Z P /Z S ) WI , from various Ward identities (WI), or else from the asymptotic value of Zf OM /Zf OM . 
We shall return to their conclusions at the end of the paper. 

However, since one has a series of values for (3, it is possible to do better; by observing the variation 
of a quantity when one increases f3, one can estimate its discretisation error separately. Indeed, the 
discretisation error will then correspond to the deviation by powers of a from what is expected close to 
the continuum : namely, as we said, one expects a cutoff independent Goldstone pole, and a very slow 
dependence of the perturbative part of Zp/Zs itself through g 2 . Note that these expectations amount, 
on the whole, to saying that Zp/Zs should be rather stable with respect to /3, on the limited range of 
available j3 values, and we will refer to this, from now on as "stability"; the discretisation errors can then 
be estimated as the deviation from this stability. 

Such a study of the (3 dependence is indeed possible for the Goldstone contribution to Zp IOM /Zg IOM 
since one has three (3's, and to some extent for Zp IOM jZ^ OM itself, although one must then take into 
account the dependence expected from 0(g 2 ) corrections, which is very slow, and one must make sure 
that the comparisons are performed for the same parameter values in physical units. 



3.1.3 Extrapolating to m q = 

In the presence of a Goldstone pole, and because of its l/m q behavior, there is no m, -» limit at all 
for l/Z P 10M , and there is the trivial one for Zf OM , or Z% 10M jZf OM . One can define a chiral limit 
only after subtracting the Goldstone pole. 

On the other hand, if, as in || one considers Zp IOM as it is, without subtraction of the pole, and 
if one then makes as usual a linear fit in m q , the extrapolation to m q = is not the chiral limit. Two 
questions then arise : 

1) Is such a linear fit possible, given that the real behaviour includes a l/m q term? 

2) What is the meaning of the quantity obtained by this linear extrapolation? 

As to question 1), a fit to l/m q linear in m q seems possible for the values of m q reached in standard 
numerical simulations, at least forp 2 not too small, but is not with smaller masses, as reached by JLQCD. 
As to question 2), we show below that in fact the linear fit used in || gives to a good approximation 
%mom i^mom at an e ff ec tive mass m e //(/3) which can be determined by the various masses used in the 
extrapolation and which is close to the lowest mass. 



3.2 Results 

First of all, the extrapolated results of can then still be used to observe the 1/p 2 power behavior of 
the Goldstone pole, and a fit in 1/p 2 will first allow us to quantify this contribution. Secondly, we shall 
determine m e //(/3) and this will enable us to observe the typical Goldstone sensitivity to the mass: the 
apparent discretisation errors are in fact due to the hidden m e //(/3) dependence. Finally, we shall be 
able to confirm the quantitative estimate of the Goldstone coefficient we have made previously, and to 
give an estimate of the true discretisation errors in Wilson data. 
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3.2.1 Large power corrections 

To display the power corrections, we perform a fit on the Wilson Z p IOM /Zg IOM ratio, given with its 

statistical errors i n M , and shown in Fig. 1. Note that the momentum variable is the true p, not 
p = sin(ap)/a as in [g|7 The values of the parameters corresponding to the data are given in Table 1. 






6.0 


6.2 


6.4 


1/a (GeV) 


2.258±0.050 


2.993±0.094 


4.149±0.161 


m 


0.1530 


0.1510 


0.1488 


K2 


0.1540 


0.1515 


0.1492 


K3 


0.1550 


0.1520 


0.1496 


K4 




0.1526 


0.1500 




0.15683 


0.15337 


0.15058 



Table 1: the parameters corresponding to the data analysed here, from 



One can clearly see in these data the presence of a pole contribution at small p 2 . To quantify it, we fit 
the points at each separately, to the form: 



TMOM IryMOM 



as/p{P) H 2~ 



(3) 



with b in GeV 2 and p in GeV. The results of this fit are shown in Fig. 1(b) and in Table 2. 

The x 2 /dof is rather high, but this is due to the points at high p 2 (this can be verified: cutting out 
high values of p 2 reduces the x 2 /dof substantially), where the discretisation errors and/or the logarithmic 
corrections should be the largest. As expected from BPT, the p 2 -independent term is remarkably stable 
with (even at 6.0), as/p — The coefficient of the power correction is less stable, changing by 

about 15%, but is consistently very large. Hence b is clearly incompatible with zero, and gives a very 
large effect of around 30% on the total zf OM /Z^ OM at 2 GeV. 






6.0 


6.2 


6.4 


a S/P 


1.1414±0.0072 


1.1266 ± 0.0082 


1.1364±0.0049 


b 


1.8710±0.063 


2.8996±0.15 


2.5844±0.15 


x'/dof 


1.49 


1.22 


1.83 



Table 2: the values of the coefficients of Eq. (g) fit to the data of [ 
The consistency of the values of b at ■ 



6.2 and = 6.4 shows that this contribution is much beyond the 
/zj P . In fact, these can be estimated to about 2% from the difference 
between Zf UM /Z^ OM at 6.2 and 6.4 at p ~ 2 GeV. 

The difference of b at — 6.2 and — 6.4 might be taken as indicating the discretisation artefact on 
the coefficient itself. However, we show in the next section that even this difference is most probably a 
physical effect, and that the real discretisation error on the power correction is still smaller. 



discretisation error on zf OM /Z MOM 

% MOM jgl 



3.2.2 Power corrections are of Goldstone origin 

Having proven the existence of large power corrections, stable with 0, and therefore probably not artefacts 
of discretisation, we must now prove that these come from a Goldstone boson. This is in agreement with 
the dominance of the divergence of axial current (the pseudoscalar density) at small pion mass, but in 
fact we can show that the data itself favours this interpretation. 

However, we would like first to comment on the dominance of the Goldstone, and, for that purpose, 
to establish the connection with the quantity Wlq discussed in ||: power corrections to Zg lOM /Zp 4014 
can originate only from the Goldstone boson pole ~ \jm q , if one is close enough to the chiral limit. 

We start with the Ward identity given in equation (19) of |0]: 



m\Y p(ap] ami, ami) — m2^p(o-P] ami, ami) 
(mi — m2)Ts{ap\ ami, ami) 



(4) 
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where Tp and T$ are the bare vertex functions. If we assume that Tp has a Goldstone contribution, 
whereas T$ doesn't, we get: 

Tp = Ap(p 2 ) + ^^ + 0(m q ) 



Ts=A s (p 2 ){l + \s(p 2 )m q ) 

We see that the r.h.s. of Eq. (g) must then be A P (p 2 )/A s (p 2 ) + 0(m q ). As (Z S /Z P ) WI is a constant, 
we have: 

with C a constant, independent of p 2 . 

But in the MOM scheme, the ratio is given by: 

Zf OM _ Y_p_ 

^mom - y; 



Hence : 



Z¥° M _ c+(m _ Mp^Bpjp 2 ) 

Z MOM - C + Ag(p 2) ^> ^2) 



We see that the power corrections are dominated by the Goldstone pole l/m q . There seems to be 
possible additional O ((m q )°) power corrections, although \s(p 2 ) is probably small^j. These are them- 
selves roughly proportional to the residue of the Goldstone term, since Xs(p 2 ) is not expected to have 
quick variation with p 2 . At least, they are connected with the presence of the Goldstone pole, and we 
can say that all the power corrections, to this order 0(m q ) included, originate in the Goldstone pole. In 
our fits, we shall neglect the Xs{p 2 ) term in Eq. (5), as well as the smaller 0(m q ) terms. 

We see also that the the difference : 

Awi - X s(p 2 )Bp(p 2 ) 

{ Zp> Z^OM- mqAsm W) +0("Vr) (5) 

is entirely due to the Goldstone boson. 
3.2.3 The effective quark mass; Goldstone fit 

If the 0(m q ) corrections are not large, the linear extrapolation to n cr i t which is usually performed 
amounts to making a linear fit in m q to l/m q , for the 3 or 4 values of m q with 2am q = 1/k— \j K cr i t . 
The extrapolation of the resulting straight line to m q = (or to K cr i t ) then defines the inverse of an 
effective mass l/m e //(y9). The extrapolated Zg IOM jZp OM is thus really calculated, not at the chiral 
limit, but at m e ff(P), and has the form: 

rrMOM y 

JMoW = a s/pW+ meff{(3)p 2> (6) 

in other words: 

(7) 



m e ff(0) 

where b' is a constant, i.e. a number independent of p, m q and /3, which we call the Goldstone strength. 

3 This (m q )° term was ommitted in the initial version of our paper. We thank D. Becirevic for having helper! ns realise 



this. That this Ag is indeed small is implied by the observations of the QCDSF group for Wilson artinn hep-lat/9807044 



p. 16 ; for Kogut-Susskind action, it is striking in the Fig 1 of the JLQCD paper, hep-lat/9901019 ; for ALPHA action, we 
thank D. Becirevic for confirming that Z$ is incredibly stable with respect to variations of m q over a very large range of 
light masses. 
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As long as the Goldstone term is not too large, this result is maintained to a good approximation 
when the linear extrapolation is made on the inverse, Zp IOM /Zg , which is what is actually done in 
§■ 

We show in Fig. 2(a) and in Table 3 the result of this extrapolation for the values of K given in Table 
1. As we see, m e //(/3) is close to the lowest mass used in the extrapolation. 



8 6.0 


6.2 


6.4 


m ef f(0) (GeV) 0.0591 


0.0400 


0.0434 


b' (GeV a ) 0.1106±0.004 


0.116±0.006 


0.112±0.006 



Table 3: the values of the effective mass defined by the linear extrapolation, and the resulting values 
of the coefficient of the Goldstone term. 
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Figure 2: (a) The three straight lines extrapolate the values of l/m q at each 0. Their intersection 
with the in = ordinate defines l/m e //(/3) and (b) the result of a joint fit to all data, after the 
l/m e ff(0) dependence of the Goldstone term has been taken into account. 



We can now deduce b' from the fit to b calculated previously: b' = 6m e //(/3). The striking result, 
shown in Table 3, is that the rather different values of b obtained previously in Table 2 correspond to 
very good approximation to the same b' , i.e., the Goldstone strengths extracted from the data are almost 
the same, and the difference of the b's is mainly due to the different values of m e //(/3) for each 8, which 
are due to the choice of values of k. 

This effect is also evident in Fig. 5 of fj and in our Fig. 1(b) , where Zf OM jZf OM at /3 = 6.0, 
with p in physical units, deviates significantly from its value at the higher /3's, at small momenta. Such 
a deviation cannot be explained by discretisation effects, which are not confined at small momenta. The 
natural explanation is that the effective mass is notably higher than at higher 8, and that the apparent 
deviation of Zp IOM /Zg at j3 = 6.0 is almost entirely due to the quark mass dependence of the 
Goldstone effect, and disappears when results are compared not only at same physical p, but also at 
same physical quark mass. The typical quark mass dependence of the Goldstone had been hidden by 
the extrapolation procedure, but has reappeared as a completely spurious discretisation effect. The true 
Goldstone origin of this fictitious discretisation effect is revealed by its l/(p 2 m q ) behaviour. 

We are now in a position to perform a joint fit to the data at the three values of 8, with the variable 
l/(m e //(j9) p 2 ) instead of l/p 2 . We find, with a x 2 /dof = 1.39: 

7 mom I7 mom n , (0.112 ± 0.025) GeV 3 

Zs /Zp ~ a s/P^) + m eff {8)p 2 (8) 
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The dependence of a s/P on (5 is very weak: we find l/a S / P ((3) = 0.88±0.05, 0.88±0.04 and 0.88±0.03 
respectively at j3 = 6.0, 6.2 and 6.4. 

The Goldstone contribution is stable, compatible at the 3 /3's within statistical errors, and very large 
when p ~ 2 GeV and m q ~ 50 MeV. Note that these results, shown in Fig. 2(b), are in direct contradiction 
with the conclusions of J(|. 

The mildness of cutoff dependence is manifest in the possibility of making such a good common 
fit to the data for the three different /3's. This possibility also gives strong support to the Goldstone 
interpretation, since it would not be possible without accounting for the l/m q dependence. 



l/ a S/p ~ 0.88, close to the 
(Ml), from which however one 



3.2.4 a 2 p 2 discretisation errors 

From the fit (||), we can deduce the asymptotic value of Zp OM /Z s AOM = 
BPT result. This should be equal to the value given by the Ward identity 
gets a lower result 0.79 - 0.80 §]. 

In fact, there seems to be a p-dependent effect, which is signaled by the fact that the strength of the 
Goldstone term becomes slightly lower and that the x 2 
Also, correspondingly, the right-hand side of Eq 



improves when cutting off the large a 2 p z points. 
[|, which equates to Zf OM /Zf OM minus the 
Goldstone, is not perfectly constant as it should, although it is much more so than Zp IOM /Z<f OM itself, 
illustrating the Goldstone interpretation (see Fig. 3 of j6|). It is also somewhat different from l/as/p, 
although the comparison between a constant and a varying quantity is difficult. 

The residual cutoff dependence (leaving aside the 0{g 2 ) effect) can be fitted by a small negative a 2 p 2 
term in Z s IOM /Z P iOM . We can then obtain very good fits at the three /3's with a universal Goldstone 
strength b' = 0.098 GeV 3 and a constant term a S /p which decreases with /3 just as expected from BPT, 
and we obtain our final result: 



gMOM i^MOM 



■ as/p{0) 



(0.098 ±0.004)Gey 3 

m e ff(P) P 2 



- (0.013 ±0.003)aV 



(9) 



with x 2 /dof = 0.45 and dof — 19. The corresponding values of a s /p, together with the expectations 
from BPT and the Ward identities, are shown in Table 3. From this table, it is visible that, once more, 
after due subtraction of the essentially non perturbative Goldstone pole effect, one has a result close to 
BPT. Note however that the BPT estimate quoted here is the ratio of the one-loop BPT estimates of Zp 
and Zs at ap = 1; one would obtain a somewhat different result, and one exactly scale independent, by 
applying one loop BPT directly to the ratio. 





6.0 


6.2 


6.4 


l/as/p ( our fit) 


0.835±0.010 


0.845±0.009 


0.845±0.007 


l/as/p (BPT) fa 


0.83 


0.84 


0.85 


1/as/P (Wlq) If 




0.79 ±0.02 


0.80 ±0.02 



Table 3: our determination of l/as/p compared with other determinations from 



We show in Fig. 3 that the Goldstone contribution dominates the discretisation artefact described by 
our a 2 p 2 term, up to the highest considered a 2 p 2 for the two higher values of /3. 
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a 2 p 2 discretisation 
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Goldstone 
contribution 

(%) 
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a 2 p 2 

Figure 3: the contribution of the a 2 p 2 term (rising curves) and of the Goldstone contribution 
(decreasing curves) relative to the total 



This fit, as shown in Table 3, also gives an estimate of the asymptotic contribution to the Ward Identity, 
1/as/p, lower and closer to (f|). We also give our estimate of the relative sizes of the three terms at 
a 2 p 2 = 1 in Table 4. 



(3 


6.0 


6.2 


6.4 


Goldstone term 


21.5% 


19.0% 


10.1% 


a 2 p 2 discretisation errors 


-0.83% 


-0.87% 


-0.96% 



Table 4: the relative values of the Goldstone boson and of the a 2 p 2 discretisation errors at a 2 p 2 — 1. 

We must of course keep in mind that the procedure used here is crude, and that a proper analysis 
of the data must take the Goldstone contribution into account for each k, before extrapolating to the 
chiral limit. Hence the agreement with BPT and the difference with the WI determination may reflect 
the crudeness of our method, due to the unavailability of better data. It is also possible to get fits still 
closer to 1/as/p = 0.8 by allowing for a 4 p 4 terms, at the cost of some variation of b 1 with (3. On the 
other hand it is clear that a similar analysis of discretisation artefacts of the type a n p n could be usefully 
applied to Wlq, which does not appear to be perfectly constant. At any rate, we observe that even with 
our latter fits with a 4 p 4 terms, the main conclusion is still the same: the Goldstone strength is large, 
close to our previous estimate, and well above any estimate of discretisation errors over a large range of 
momenta above 2 GeV. 

3.2.5 Comparison with Goldstone residue extracted from QCDSF improved data 

We can compare the Goldstone residue with our result of Eq. (||), extracted from the QCDSF data. 
Since we do not know the value of Zs from the ALPHA action, the simplest thing to do, disregarding 
slow logarithmic evolutions, is to compare the magnitudes of the ratio of the power correction to the 
perturbative term at p and m q similar in physical units, or of the Goldstone strength b' to the perturbative 
term in the same unit GeV 3 . We immediately notice that the latter ratio is in perfect agreement: 0.82 
GeV 3 for the improved action, 0.818 ± 0.003 GeV 3 for the Wilson action. 

Hence there is full compatibility between the various determinations of the Goldstone strength, both 
at various /3's and for various actions, converging towards very large values, some 30% of the total around 
2 GeV and for a mass around 50 MeV. 
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4 Discretisation errors; an overall discussion 



We now come to the paper [j| of Giusti and Vladikas (G&V), and to a discussion of the origin of 
its conclusions, opposite to ours. Admittedly, many of the procedures used there are common in the 
literature, but they turn out to be inappropriate for the present discussion of the Goldstone contribution. 
Hence we think that beyond answering the criticisms of commenting upon them is of general interest. 

Although they intended to discuss specifically our estimate of the magnitude of the Goldstone contri- 
bution, the core of the argument of G&V is the comparison of various determinations of (Zp/Zs) WI . The 
spread of the values directly extracted from Ward identities, and the deviation with the values obtained 
for Z^ OM /Zf OM , are supposed to be a measure of the discretisation error on Zf OM /Zf OM itself. 
This way of estimating the errors is one of our main disagreements, for reasons expressed below in points 
2 and 3. The other main difference is that G&V do not take into account the strong scale dependence 
of Zf OM /Zf OM ', as explained in points 1 and 4. 

1. Let us first emphasize that, even admitting the 10 — 15 % estimate of discretisation errors made by 
|Q, these errors cannot dominate the 30 % estimate of the Goldstone pole that we have given at 
2 GeV. The reason why G&V have missed this point is clear. The argument, illustrated in their 
Fig. 4 and Table 1, relies on a point with large momentum for each (3. The numbers of ref. || are 
given at sin 2 (ap) = 0.8, i.e. at p = 3.3 GeV (/3 = 6.2) or p = 4.6 GeV (/3 = 6.4). This choice of 
a large physical momentum is not appropriate when the manifest goal is to discuss the Goldstone 
pole overall strength. Indeed, if there is a Goldstone pole, Zp OM /Zg IOM strongly depends on p 2 , 
and its difference with the asymptotic value decreases rapidly with increasing p 2 , rendering difficult 
or eventually impossible the determination of the power correction. Had G&V taken p = 2 GeV, 
they would have had to quote a central value for Zp IOM /Zg IOM around 0.6 or less (as can be seen 
from their Fig. 3), much below Wlq = 0.79 and also below Wlh — 0.68 — making manifest the 
large magnitude of the Goldstone — . In the introduction to their new version of the paper ||, they 
state however that Zp IOM /Zg IOM is compatible with the WI's even at 2 GeV within discretisation 
errors; this, from their own numbers, amounts to admitting still larger discretisation errors of the 
order of 30 % or more at j3 = 6.2 and a 2 p 2 — 0.45 (difference between 0.79 and 0.6). One the 
other hand, the only known way to explain the data with a reasonable error estimate is through 
the Goldstone interpretation. 

2. Furthermore, let us emphasize that their estimated discretisation errors contradict the evolu- 
tion of data with j3 as far as Z P IOM /Zg IOM is concerned. Indeed in [Q, the statistical errors 
on Zp OM /Zg IOM are small, and the discretisation errors seem also small, since the values of 
Zp IOM /Zg IOM taken at the same physical momenta differ only by a few percent between (3 = 6.2 
and (3 = 6.4 (see Fig. 1(a)). 

A careful reading of the text reveals that the much larger error introduced in J(| has actually nothing 
to do with the error on Zp IOM /Zg IOM itself, but really concerns the estimated error made on the 
indirect estimate of the Ward identity result through the asymptotic value of Zp IOM /Zj' fOM . This 
error was already given in jjj, and in fact, in 0, the same numbers were quoted as "RGI" (i.e., 
estimate of the asymptotic, renormalisation group invariant quantity). Indeed, in jj, the lack of 
the expected plateau was interpreted as an error of 10 to 15% on (Z P /Z S ) RGI . Intel, the same 
number is re-expressed arbitrarily as an error on the value of Z^ OM /Z s AOM at the lower end of 
the range, a 2 /j, 2 = 0.8, though it is not an actual error on Zp IOM /Zg IOM . 

In our opinion, the lack of plateau signals power corrections and the need to subtract them. The 
procedure of ||, which amounts to including them automatically into the errors on Zp OM jZ s aOM , 
makes it of course impossible to discuss the power corrections. 

3. Moreover, G&V further substantiate their estimate of the discretisation effect by observing a 10 — 
15% discrepancy between determinations from two Ward identities, called Wlq and Wlh, at (3 = 6.2. 
This procedure has the advantage that the Ward identities are scale independent. It is also a natural 
approach, if one is working at only one (3, to look for the difference between quantities which should 
be equal. However, in this approach, one does not know which is the best estimate, or whether 
both equally fail: the same discretisation error is attributed to both, and to any other quantity, 
such as (Zp/Zs) RI / MOM , which may be over-pessimistic. 
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As already emphasized, a better approach is to examine the variation of the specific quantity 
which one wants to study, when one increases (3. Admittedly, small variations can present from 
the variation of as, but they should be 0(g 2 ) in BPT and vary slowly with (3, and hence the 
quantity should be very stable when (3 changes. If it is stable up to logarithms, this particular 
quantity has probably a small discretisation error. It seems that the Wlq determination changes 
very slightly between 6.2 and 6.4, from 0.79 to 0.8. In fact, the values are compatible within 
statistical errors, and the small increase is expected from BPT. This is not so for Wlh, which shows 
a strong variation from 0.68 to 0.73. In fact, Wlh corresponds to Za times the ratio p/m q of 
the axial to the subtracted mass, and the latter ratio is known to exhibit rather large variations. 
The natural conclusion would be then that Wlq deserves more trust than Wlh. The same can be 
said probably of Z^ OM /Zg IOM , which is remarkably stable, with a small variation in agreement 
with BPT. . Thus the large discretisation error should be probably attributed to Wlh only, not 
to the three quantities at the same time. It is then rewarding that Wlq and Z P IOM jZ^- OM give 
compatible results for the estimate of WI, after due subtraction of the Goldstone pole, as we have 
shown above. 

4. We note that in their Fig. 4 and Table 1, G&V consider the spread of values of various determi- 
nations of Z P /Z S with increasing (3, including Zp IOM /Zg IOM , at the same a 2 ^ 2 = 0.8, and not at 
at the same physical p, as one should do when discussing the error on a scale-dependent quantity. 
G&V are in fact not comparing the same quantity at two different f3's, but two different quantities: 
the values of zf OM jZf OM respectively at p = 3.3 GeV and p = 4.6 GeV. 

The natural explanation of the increase of Zp IOM /Zg IOM with f3 in their figure and table is the 
decrease of the power correction with increasing p, which is a physical effect, not the discretisa- 
tion errors, except at very large p. If we duly compare at identical p, we see once again that 
Zp IOM /Zg tOM is very stable with (3, and that the discrepancy with Wlq does not decrease, unlike 
suggested by the Fig. 4 of ||: it is a physical effect, the sign of Goldstone contribution as shown 
above. 



5 Conclusion 

The quantity Z P IOM /Z^ OM appears to be the best indicator of the Goldstone pole. Contrarily to ||, 
we find a large Goldstone contribution to the Wilson data, of the same magnitude as found previously 
with data for Zp IOM from the QCDSF improved action. The results are consistent for 3 values of (3 
and for momenta ranging from about 1 GeV to 8 GeV. Of course, the determination of the Goldstone 
strength comes mainly from moderate momenta, where the contribution is the largest. The discretisation 
uncertainty, as estimated from the variation of Zp IOM /Zg IOM with (3, appears in fact to be rather small, 
and the Goldstone contribution at p — 2 — 4 GeV is far above it. Even admitting the larger discretisation 
error advocated by which is not relevant in our opinion, our claimed Goldstone contribution at 2 GeV 
is so large, as already found previously, that it is clearly dominating. Evidently, it is smaller at the higher 
momenta considered by G&V in their Fig. 4 and Table 1, but this is as it should be: it must be oc 1/p 2 ! 

This large Goldstone contribution explains in a natural manner the discrepancy of the MOM Zp with 
the ALPHA group determination of Zp at large distance (around 30% at 2 GeV). It also explains for 
the most part the absence of plateau in Z p 10M /Zg I0M , even at the highest acceptable momenta. True, 
we find some contribution from a 2 p 2 artefacts, but certainly not a dominant one. Given this absence of 
a plateau, one should not insist on extracting an estimate of (Zp/Zs) WI directly from Zp OM jZ^ OM , 
even if one assumes large errors. The only way to proceed, which we have illustrated here, is to subtract 
the Goldstone contribution. The result of the subtraction is, once more, a number close to the BPT ex- 
pectation, which is quite encouraging. Another formulation of this is to use Eq. (Q), which automatically 
subtracts the Goldstone, and which is found to be rather stable with {3. 

Of course, some slight changes in the conclusions must be expected from a more thorough analysis 
of the complete data, where it may be possible, in particular, to explain the small discrepancy between 
Wlq and \ J ag jp given in Table 3. 

It remains to explain the apparently different conclusion from [jn^, at 13 = 6.2, which find smaller 
power corrections; this may be related to off-shell improvement. 

The interesting and intriguing physical question is now to find the reason why the Goldstone residue is 
so large in Zf OM or Zf OM jZf OM . This is connected with the behavior of the pion wave function (BS 
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amplitude) at short distance. The apparent contradiction with the standard OPE is puzzling. A naive 
interpretation of our results^] would be to claim that the quark condensate is 10 times the standard value, 
but this is certainly not probable, and the ultimate physical reason of this disagreement must surely be 
more subtle. 
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Note added in proof: 

While we were writing this letter, C. Dawson [Q and Y. Zhestkov ju|] stressed again the necessity of 
subtracting the Goldstone pole to obtain a chiral limit. 
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